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Abstract 

Ab initio gap equation for 1 Sq pairing in a nuclear slab is solved 
for the Argonne vis iViV-potential. The gap function is compared 
in detail with the one found previously for the separable form of the 
Paris potential. The difference between the two gaps turned out to 
be about 10%. Dependence of the gap on the chemical potential fi is 
analyzed. 
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1 Introduction 



In the last decade, a progress was reached in the microscopic theory of pairing 
in atomic nuclei. In this connection, the studies of the Milan group [H [2] 
should be cited first. In the recent paper [2] the direct solution of the ab 
initio gap equation for a nucleus 120 Sn was carried out using the realistic 
Argonne Vu iViV-potential. A bit later we made a similar calculation [3] for 
a slab of nuclear matter with parameters mimicking the atomic nuclei of the 
tin region with the Paris potential pi]. The comparison of results of these two 
calculations rises some questions. In Ref. [2] the gap A, obtained by solving 
the gap equation, turned out to be about a half of the experimental value. 
The discrepancy was attributed to the missing contribution of the low-lying 
surface vibrations [2]. On the other hand, the gap found in [3] is much closer 
to the experimental one and leaves a room to the surface vibration corrections 
only for a 20% about. 

To some extent, the difference could originate from the difference between 
the two systems under consideration, the spherical nucleus and the slab. This 
reason seems to us not very important as far as only the curvature effects 
make the difference. Different values of the effective nucleon mass m* used 
in these two alternative calculations look more essential. Indeed, it is well 
known [5] that the gap value in infinite matter is rather sensitive to the 
value of the effective mass m* . In [2] the coordinate dependent effective 
mass m*(r) was used of the Sly4 version [6] of the Skyrme force. Its value 
changes from m* = 0.7m inside a nucleus to m* = m outside. In [3] the 
bare nucleon mass, m* = m, was used, in accordance with prescriptions of 
[7] and [8]. The specific form of the realistic AW-potential used in the gap 
equation is one more possible reason of the difference under discussion. In 
this article we analyze the latter point repeating the calculation of [3], but 
for the Argonne Vis potential [9] which is quite similar to that of [2] but 
essentially different from the Paris potential. Indeed, both the version of 
the Argonne force are rather soft core whereas the Paris potential is one of 
the most hard core realistic AW-potentials. In addition, we examine the 
dependence of the pairing gap on the chemical potential fi of the system 
under consideration. 



2 



2 The formalism 



The general method for solving the gap equation in a system with the ID 
inhomogeneity is developed in Ref. [3]. Here we present only the formalism 
which is necessary to describe the procedure. Just as in Ref [3], we use the 
two-step method [10] to solve the microscopic gap equation, 

A = VA S A, (1) 

where V is the free iViV-potential and A s = GG S stands for the two-particle 
propagator in the superfluid system. Here G and G s are the one-particle 
Green functions without and with pairing effects taken into account, respec- 
tively. In Eq. (pQ), as usual, integration over intermediate coordinates and 
summation over spin variables is understood. In the nuclear matter theory, 
the gap equation in the form of Eq. ([1]) is usually referred to as the Bardeen- 
Cooper-Schrieffer (BCS) approximation^ In the more profound many-body 
theory [TTI [T2] the irreducible block of the iViV-interaction should replace V 
in Eq. <^). 

The complete Hilbert space S of two-particle states is split into two parts, 
S = So + S'. The first one is the model subspace So in which the gap 
equation is considered, and the other is the complementary subspace S'. 
They are separated by the energy E in such a way that S involves all 
the two-particle states (A, A') with the single-particle energies E\,E\r < E . 
The complementary subspace S' involves such two-particle states for which 
one of the energies s\,e\/ or both of them are greater than Eq. Therefore, 
pairing effects can be neglected in S' if E is sufficiently large. Validity of 
the inequality A 2 /(E — fi) 2 <C 1 is the criterium of such approximation. 

Correspondingly, the two-particle propagator is represented as the sum 
A s = Aq + A'. Here we already neglected the superfluid effects in the S'- 
subspace and omitted the superscript "s" in the second term. The gap equa- 
tion ([1]) can be rewritten in the model subspace, 

A = ^A, (2) 

where the effective pairing interaction (EPI) should be found in the supple- 
mentary subspace, 

v: s =v+va'v: s . (3) 

In physics of finite nuclei this term is commonly used for an approximate way to solve 
the pairing problem, contrary to solving the Bogolyubov equations. 
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As is well known, the gap equation ([T]) can be written in an equivalent 
form, 

A = -Vx, (4) 

where the anomalous density matrix x is expressed directly in terms of the 
Bogolyubov u-, v-f unctions: 



x(ri,r 2 ) 



J2 u i( r i) v i( r 2)- ( 5 ) 



The summation in (jSJ) is carried out over the complete set of the Bogolyubov 
functions with eigenenergies Ei > 0. 

The gap equation (TS]) in the model subspace can be also written in a form 
similar to 

A = -KVo , (6) 
where x is the anomalous density matrix specified in the model subspace as 

x (ri, r 2 ) = J A S (E; r h r 2 , r 3 , r 4 )A(r 3 , r 4 )dr 3 dr A . (7) 

However, a simple expression for this quantity similar to Eq. ((Sj) does not 
exist. To obtain the explicit form of xo it is convenient to use the basis of 
the single-particle functions <fi\ (r) which makes the normal Green function 
G diagonal: G\y(e)=G\(e)5\\> . Then one gets: 

u t (r) = J2^Mr), (8) 

A 

^(r) = 5>, A A (r), (9) 

A 

A(r 1; r 2 ) = A A y<Mri)<Mr 2 ). (10) 

A, A' 

After a simple algebra [3], one finds 

(o) 

x ( ri ,r 2 )= K^^ + a-nAj^ 2 ^ 1 )^^)^^), (11) 

i,Ai,A2 

where rt\ = (0; 1) are the occupation numbers without pairing. The super- 
script "(0)" in the sum means that all the states A belong to the model 
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subspace: S\ < E . It can be easily seen that in the limit of E — * oo, when 
the set {A} is complete, the expression (TTTT) coincides with Eq. (j5J). 

We consider a slab of nuclear matter embedded into the one-dimensional 
potential U(x). In this case, the momentum kj_ in the (y, ^)-plane (or s- 
plane) is an integral of motion. The single-particle wave functions can be 
written as 

A (r) = e^*y n (x), (12) 
and the Bogolyubov functions have the form: 

Ui(r) = e l ^ s u t (k ±J x) (13) 

v i (r) = e ik ^v i (k ± ,x), (14) 

Just as in Ref. [3], we will use the mixed representation, i.e., the co- 
ordinate representation for the x-direction and the momentum one for the 
s-plane. The scheme of solving the gap equation in the mixed representation 
was developed in detail in [3] for the case of the separable form of the Paris 
potential [131 E]- Now we consider the original form of the Argonne vi§ 
potential [9] and present the main formulas for the general case. 

Let us begin with the EPI. In the Wigner representation it depends on 
the relative momenta ki and k2 and on the CM coordinates X\ and X 2 in 
the input and the output channel. As far as we deal with the singlet pairing, 
the zeroth harmonic of the EPI enters into the gap equation which depends 
only on the absolute value of the relative momenta. One more advantage of 
the two-step method is the possibility of using the so-called local potential 
approximation (LPA) [10] to find the EPI. According to the LPA, the EPI 
at each point X = (X\ + X 2 )/2 can be replaced by the one in infinite nuclear 
matter placed in the external potential field U(X): 

V™ A (h, k 2 ; X h X 2 ) = V^(h, k 2 ; t, [U(X)\), (15) 

where t=X\ — X 2 . To find the EPI in infinite system, it is convenient first 
to solve the Lippman-Schwinger equation in the complementary subspace in 
momentum space, 

>S(M„P;%v( t „ fe ) + / (16) 

JkeS' (^r E - £( p > k ) - £ ( p > k ) 

where E = 2fi, e(P, k) = (P/2 + k) 2 /(2m) + U(X), e'(P,k) = (P/2 - 
k) 2 /(2m) + U(X). The integration in ffTBT) is carried out over k G S'. The 
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quantity V'g (ki, k 2 ; t, [U(X)\) is found from V^{k u k 2 ; P, [U{X)\) with the 
inverse Fourier transformation. 

The explicit form of the gap equation ([6]) in the Wigner representation 
for the ID geometry is as follows: 

A(k,X) = - j V^(k,k f ;X,X') Mk' ± ,k' x ;X>) ^jL^kdX' , (17) 
The anomalous density can be written as 

(0) 

Mk±,h;X) = J2< n '(k ± ) f nn/ {k x ,X) , (18) 

nn' 

where 

fnn'(k x ,X) = J e- lk * x y n (X + x/2)y n ,(X-x/2)dx, (19) 

and 

(o) 

< n '(k ± ) = J2 (n k±n vt(k ± )v?(k ± ) + (l-n hxn )vf(k ± )v?(k ± )) . (20) 

i 

Here we used the expansion of the Bogolyubov functions into the basis 

{VnY 

(0) (0) 

Ui(k±, x) = y n (x) , Vj(k±, x) = y~]vj{kx) y n {x) . (21) 

n n 

In this representation, the Bogolyubov equations have the form: 
{e nkx - n) u?(k ± ) + Y, A nn'(k ± ) vf(k x ) = E lU Uk±) , 

n' 

J2&nn>(k ± ) uf(k ± ) - (e nk± - n) v?(k ± ) = E iV «{k x ) , 

n> 

where e n k ± =e n + k± 2 /2m and 

A nn ,{k ± ) = J A(k ± ,x,X)y n (X + ~)y n ,(X-^)dxdX = 

/dk 
A(k,X) f nn ,(k x ,X)dX — (23) 
2n 



(22) 



The chemical potential /x in (1221) is determined inverting the equation for 
the number of particles per a unit surface of the slab: 



The above relations determine the gap function in the slab system. 

3 Calculation results 

Just as in [3], we use the Saxon- Woods potential U(x) symmetrical with 
respect to the point x = with typical for finite nuclei values of the potential 
well depth Uq = —50 MeV and diffuseness parameter of d — 0.65 fm. For the 
main part of calculations the thickness parameter of the slab was chosen as 
L = 6 fm to mimic nuclei of the tin region. As far as the Hamiltonian of the 
system is symmetrical under the axis reflection x —>■ —x, the eigenfunctions 
y n can be separated into even, ?/+, and odd, y~, functions. It simplifies 
the numerical procedure essentially, as far the n, n' indices in (|22[) should 
numerate the states of a fixed parity. 

The calculation scheme is rather different from that used in [3] for the 
separable form of the iViV-force. The common feature is the convenience 
to use the scheme of parallel calculations for different k± momentum, since 
the slab is homogeneous in the s— plane. The total number of processors we 
used is 210, which gives required accuracy for calculations with the model 
space 5*o specified by E = 15 MeV. The discrete spectrum representation 
method was used to take into account the continuum, < E\ < Eq. The 
infinite wall was put at the distance R = 25 fm. The independence of 
results on the value of R was checked. The gap equation (fT7|) was solved 
with iterations. The procedure was stopped when the maximal difference 
of values of A(k,X) obtained in two subsequent iterations became less than 
10~ 6 MeV. The typical number of effective iterations was about 10. The 
chemical potential was found anew according (I24p on each iteration. The 
self-consistent value of the chemical potential /i= — 7.95 MeV is very close 
to the one, /io= — 8 MeV, without pairing. The difference SfM— — 0.05 MeV 
is in agreement with the standard estimate 5 fx ~ —A 2 /e-p- 

To illustrate the solution graphically, we calculated the "Fermi averaged" 




(24) 



gap 



A F (X)=A(k F {X),X). 



(25) 
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Here the local Fermi momentum is defined as follows: kp(X)—2m(fi— U(X)) 
if 2m([i—U(X)) > 0, and A;p(X)=0 otherwise. Analogously, we could define 
the Fermi averaged EPI as 

V e F ff (X) = J dtV: s (h=k 2 =k F (X);t,X). (26) 

There are two parameters which influence essentially difficulties of the 
calculation procedure. They are the separation energy Eq which defines the 
model space Sq and the cut-off momentum K max in the integral of Eq. (TIB"]) . 
Let us begin from choosing the first one. Fig. [I] shows the dependence of 
the solution on the value of Eq. Analyzing these curves one can see that the 
variation of E from 5 MeV to 10 MeV results in the change of the maximum 
of the Fermi-averaged gap function of the order of 5% whereas the next step 
to 15 MeV diminishes this difference to 1%. A more quantitative information 
could be obtained from the matrix elements of the gap. The diagonal matrix 
elements A nn (k±=0) for different E are presented in Table 1. The even 
states are located in the upper part of the table, the odd ones, in the lower 
part. One can see that the convergence of A nn values increasing Eq is slower 
than that of the maximum of A-p(X). Now the average change of a matrix 
element with the variation of Eq from 10 MeV to 15 MeV is about 3%. This 
is, of course, the upper estimate of the accuracy if we choose Eq = 15 MeV 
as an appropriate value. Analyzing the convergence of A nn values in Table 
1, one can suppose that the change under consideration will be about 1-2% 
at the next step to Eq = 20 MeV. Therefore more realistic estimate of the 
accuracy of the calculations for Eq = 15 MeV is 1-2%. 

The latter could be referred to both approximations used in the version 
of the the two-step approach outlined in Sect. 2. These are: neglecting of the 
pairing in the complementary space S' and using the LPA as well. Note that 
in the case of the Paris potential 1% accuracy was reached at Eq = 20 MeV 
®. 

Let us turn now to the analysis of dependence of the gap function on 
the cut-off momenta K max in Eq. ( TTB"j) . The Fermi-averaged gap A F (X) is 
displayed in the left panel of Fig. [2] for different K max in the case of the 
Argonne Vi$ potential. The value of K max = 6.2 fm" 1 corresponds to the 
cut-off energy E max = 800 MeV used in [2]. We see that this value of K max 
guarantees only 10% accuracy of the gap value and only K max = 10 fm -1 is 
sufficiently big to result in 1% accuracy. Therefore it is possible that the gap 
found in [2] is underestimated at approximately 10% due to not sufficiently 
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Table 1: Diagonal matrix elements A nn (k ± =0) (MeV) in the slab with the 
width parameter L=Q fm for Argonne v 18 potential for different Eq (MeV). 
e n (MeV) is the single particle spectrum. 





A nn 




E = 


E = 5 


E = 10 


E = 15 


-48.26 


1.516 


1.362 


1.339 


1.284 


-36.55 


1.461 


1.282 


1.245 


1.210 


-18.60 


1.198 


1.034 


0.982 


0.972 


-1.325 


0.421 


0.342 


0.324 


0.309 


-43.51 


1.522 


1.339 


1.309 


1.266 


-28.04 


1.353 


1.174 


1.119 


1.102 


-9.10 


0.903 


0.764 


0.747 


0.717 



big value of -E max - Note that another version of the Argonne force, vu, 
was used there but it is quite close to that used by us. In the systematic 
calculations, the values of Eq = 15 MeV and K max = 10 fm -1 will be chosen. 

For comparison, the analogous dependence is shown for the Paris potential 
in the right panel of Fig. [2] demonstrating very slow convergence of the inte- 
gral in Eq. ffTBT) in this case. Indeed, now the huge value of K max = 160 fm -1 
is necessary to guarantee 1% accuracy for the Paris potential. Only the use 
of the separable version of the Paris force in [3] made it possible to carry out 
calculations with sufficient accuracy. 

Fig. [3] demonstrates that the EPI found from Eq. ( {TBI depends only a 
little on the specific form of the realistic iViV-potential used provided the 
convergence of the integral in the momentum space is reached. One can see 
that the two curves for different potentials are nearly coincident. To show 
a small difference we displayed the same curves at small and large values 
of x separately in a large scale. We see that the difference is about 2% in 
the asymptotic region where the EPI is big and about 5% where the EPI 
is small. However, even such a small variation of the EPI is not negligible 
as far as any change of the pairing force is significantly enhanced in the gap 
equation. As the analysis of [3J has shown, in the system under consideration 
a 1% variation of the EPI results in approximately 5% variation of the gap. 
It is confirmed with direct comparison of the gap functions for both the 
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Figure 1: The Fermi averaged gap Ap(X) for different Eq. 

potentials which is presented in Fig. HI The maximual difference between 
values of Ap(V) for different potentials is of the order of 10%. 

To analyze this difference quantitatively it is instructive to compare the 
matrix elements of the gap. They are presented in Table 2 for the diagonal 
matrix elements A nn (k±=0). One can see that the difference under consid- 
eration is again of the order of 10%. The only exception is the state with 
very small energy e n = — 1.3 MeV for which the matrix element of A itself is 
essentially less than the typical values which is of the order of 1 MeV. This 
peculiarity of states with small energy is caused with a special form of their 
wave functions which have very long "tails" outside the slab. Therefore the 
weight of y n {x) 2 in the integral of A nn is rather small in the region where 
the A(X) is large. The upper half of the table contains the positive parity 
states. These are the matrix elements which could be related to those in 
heavy nuclei over the single-particle s-states. Indeed, the k± variable in a 
slab is an analogue of the orbital angular momentum I in a spherical nucleus, 
k± = corresponding to / = 0. The average value of A nn (k±=0) for the 
positive parity states is A=0.94 MeV for the Argonne force and =1.04 MeV 
for the Paris potential. It looks reasonable to exclude the state with anoma- 



10 




Figure 2: The Fermi averaged gap A-p(X) for different cut-off momenta 
K max (in fm _1 ) for the Argonne Vis potential (the left panel, Eq = 15 MeV) 
and the Paris potential (the right panel, E = 20 MeV). 



lously small energy for calculating the average value of the gap. In this case, 
one finds A=1.16 MeV for the Argonne force and =1.27 MeV for the Paris 
potential. These values should be compared with the experimental values 
of the gap in the tin region which are of the order of A cxp ~ 1.3 MeV. We 
see that although the gap value for the Argonne force is a little (10%) less 
than that for the Paris force it is also rather close to the experimental data, 
reserving only narrow room (~ 10 — 20%) for contributions of the surface 
vibrations. 

The quantity Ap(X) is a localized representation of the gap. To ex- 
amine the nonlocal properties of the gap, it is reasonable to calculate the 
X-dependent rms radius 



r*(X) = yViPO), (27) 



where 



(ri(X)) = 



q2 

Jr 2 A(r,X)dh _ -Q^ A ( k > X ) 



jA(r,X)d 3 r ~ A(k,X) 



(28) 

k=0 
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Figure 3: The Fermi averaged EPI V^ S (X) for the Argonne 
Vis potential (solid line) and the Paris potential (dashes), 
E = 20 MeV) for both the cases. 



It is displayed in Fig. [51 We see that it is about 0.8 fm inside the slab 
and grows to ~ 1.2 fm outside, i.e. the gap function is essentially non-local. 



4 On the /i-dependence of the gap 

In this section, we analyze dependence of the pairing gap on the chemical 
potential /i. This subject is of interest due to a progress in physics of of nuclei 
laying far off the /^-stability valley reached in recent years. In particular, it 
concerns nuclei in the vicinity of the nucleon drip line, where the chemical 
potential \i of the neutron or the proton subsystem vanishes. The importance 
of nucleon pairing for this region is generally accepted. In studies aimed at 
calculating the nucleon drip line within the Hartree-Fock method [151 EE] or 
the generalized EDF method [5] , which employ a phenomenological nucleon- 
nucleon interaction, much attention was given to exploring special features 
of pairing at small values of /i. In this connection, the main interest was 
put to such problems as the correct taking into account continuum states, 
as a comparison of the exact solution of the Bogolyubov equations with that 
found within the BCS approximation, and so on. However, the possibility 
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Table 2: Diagonal matrix elements A nn (k ± =0) (MeV) in the slab with the 
width parameter L=6 fm. e n (MeV) is the single particle spectrum. 









Argonne Vi 8 potential 


Paris potential 


-48.26 


1.284 


1.396 


-36.55 


1.210 


1.326 


-18.60 


0.972 


1.077 


-1.325 


0.309 


0.357 


-43.51 


1.266 


1.385 


-28.04 


1.119 


1.227 


-9.10 


0.764 


0.817 



t 1 1 1 1 1 r 




X,fm 

Figure 4: The Fermi averaged gap A-p(X) for the Argonne V\s potential 
(solid line) and the Paris potential (dashes). 

that the parameters of the EPI determining A change in the vicinity of the 
drip line has been ignored so far. It is hardly possible to address this issue 
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Figure 5: The rms radius of the gap A-p(X) for the Argonne vis potential. 

within a phenomenological approach. This requires the use of the ab initio 
methods which start from a free NN-interaction. 

Such a study was carried out in jTTJ for the Paris potential. Dependence 
of A on n turned out to be model dependent provided the mean field U(x) 
is not found self-consistently. The matter is that, as it was shown in p^8l [19] . 
the potential well U(x) could change in vicinity of the drip point significantly 
due to dependence of the scalar Landau-Migdal amplitudes [11] on \i. As far 
as this problem is not solved up to now, two models were examined in (TTj. 
In the first one (Model 1), the depth of the potential well Uq=— 50 MeV is 
fixed, the central density increasing with decrease of In the second model 
(Model 2), the central density is considered to be a constant, the value of 
\Uq\ decreasing with decrease of \fi\. It turned out that, with decrease of 
the gap value decreases in the first model and increases in the second one. In 
addition, the \i dependence of the so-called gap-shape function was examined 
in [T7|. This quantity is defined as the ratio 



Xf(X) 



A F (0) • 



(29) 
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Figure 6: The Fermi averaged gap Ap(X) for different values of the chemical 
potential \i. 

It determines directly the effect of the surface enhancement of the gap [2U1I2T]. 
It turned out that this effect grows with decrease of within both the 
models. 

Now we carry out analogous calculations for the Argonne force vig. De- 
pendence of the gap on the chemical potential for both the models is displayed 
in Fig. [6j In general, it repeats the behavior of the gap in the case of the 
Paris potential [17]. The main origin of the effect is the /i-dependence of 
the EPI which is shown in Fig. [7J In the first model, with decrease of 
the EPI becomes a little stronger in the asymptotic region and weaker inside 
the slab. The first effect originates from closeness to the pole in the free 
iViV-scattering amplitude, the second one, due to growing of the local value 
of the Fermi momentum. It turned out that, with diminishing of the 
second reason dominates and A falls, and only at very small |/i| the surface 
enhancement of the EPI overcomes. In the second model, the EPI inside the 
slab doesn't practically depend on /i, and only the surface effect remains. 
Therefore in this case the gap value is growing with decreasing, but the 
effect is rather small. 

In the first model, the /i-dependence effect is much more pronounced 
for the gap-shape function, see Fig. [HJ left panel. In this case, the surface 
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Figure 8: The gap shape function xPO f° r different values of the chemical 
potential /i. 



enhancement of the gap becomes stronger with decrease of the effect 
being monotonic. It occurs because the nominator in ([2"9l diminishes slower 
than the denominator. In the second model, the denominator in (1291) al- 
most doesn't depend on /i, and qualitatively the /^-dependence of gap-shape 
function repeats that of the gap itself. 
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5 Conclusions 



We solved the microscopic gap equation for 1 Sq pairing in a nuclear slab for 
the Argonne vi 8 iViV-potential. The slab is embedded into the Saxon- Woods 
potential well U(x) with parameters typical for heavy nuclei. The width 
parameter L = 6 fm was chosen to mimic nuclei of the tin region examined 
in [2] . The features of the gap function obtained are very close to those found 
previously for the separable form of the Paris potential [3] , but the absolute 
value of A turned out to be 10% less. This difference should be attributed to 
essentially different form of the free iVW-potentials used which is the main 
input of the microscopic gap equation. Evidently, it may be considered as 
the upper estimate of the uncertainty originating from the specific choice 
of the iViV-potential. Indeed, the two potentials are very contrast to each 
other, the Argonne force being rather soft and the Paris one, very hard. The 
diagonal matrix elements of A obtained are about 1.1 MeV which could be 
compared with the experimental value of A ~ 1.3 MeV in tin isotopes with 
A ~ 120. On the other hand, the gap value obtained is significantly bigger 
than A ~ 0.7 MeV found in [2] where the lack of the gap was attributed to 
the contribution of the surface vibrations. Omitting the discussion of possible 
uncertainties in calculating these corrections (see, e.g., [10]), we mention only 
that the contribution of about 20% was obtained in [22] for the process under 
consideration. This estimate agrees with our calculations. 

We think that different values of the effective nucleon mass m* used is 
the main reason of the contradiction under discussion. In [2] the coordinate 
dependent effective mass m*(r) was used from the Sly4 force [6] which is 
essentially less than the bare nucleon mass, m* = m, used by us. The latter 
is chosen in accordance with prescriptions of [7] and [8] which are based, 
in particular, on the analysis of the single-particle spectra of magic nuclei. 
Indeed, it is well known that the use of m* < m in any Skyrme-type forces 
results in the single-particle spectra which are significantly decompressed. 

We analyzed also the dependence of the gap on the chemical potential 
/i. The result turned out to be model dependent in the approach with the 
potential well U(x) considered as an input of the calculation scheme. In the 
first model, the depth Uq=— 50 MeV is fixed, the central density increasing 
with decrease of \fi\. In this case, the value of A decreases with diminish- 
ing of |/x| and only at very vicinity of the drip point it begins to grow. A 
different behavior was found of the gap-shape function which characterizes 
the surface enhancement of the gap. This quantity turned out to be mono- 
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tonically increasing with decrease of \fi\. In the second model, the central 
density is supposed to be constant, the depth of U(x) decreasing with de- 
crease of \fi\. In this case, the /^-dependence effect is less pronounced but the 
gap value and gap-shape function grow both with decrease of \fi\. Thus, a 
self-consistent calculation of the potential well depending on /i is necessary to 
clarify the properties of nuclear pairing in the drip line vicinity. In this case, 
a /i-dependence of the effective force should be taken into account [TS], [TS] if 
one deals with the phenomenological approach. 
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